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Abstract. - The thermal conductivity k of the spin-1/2 chain cuprate BaCu2Si207 was 
measured along different crystallographic directions in the temperature region between 0.5 and 
300 K. The thermal conductivity along the chain direction considerably exceeds that along per- 
pendicular directions. Near the antiferromagnetic transition at Tjv = 9.2 K the data indicates 
enhanced scattering of phonons by critical fluctuations in the spin system. A comparison of 
the data above Tjv with available results on similar materials reveals similarities in the main 
features of the temperature dependence of the mean free path of itinerant spin excitations. 
This universal behavior is most likely caused by the spin-lattice interaction. 



According to theoretical predictions, the energy transport in several types of S* = 1/2 
one-dimensional (ID) antiferromagnetic (AFM) spin systems is expected to be remarkably 
different from similar three-dimensional (3D) systems. In 3D systems, the thermal conduc- 
tivity of magnetic origin can be large only much below the long-range ordering transitions, 
where the itinerant spin excitations are well-defined quasiparticles (magnons). At tempera- 
tures close to but below the 3D ordering transition, the spin-mediated thermal conductivity 
Kg is considerably reduced because of enhanced magnon-magnon scattering 012]. Above the 
transition, the energy transport by spin excitations is small and has diffusive character, such 
that in real materials other carriers of energy, most notably phonons and, in electrical conduc- 
tors, electronic quasiparticles, dominate the heat transport. In strictly ID spin systems, no 
long-range order exists at nonzero temperatures, and naively one could expect a rather small 
thermal conduction via the spin system. However, for ID = 1/2 spin systems described by 
the Heisenberg AFM XX Z Hamiltonian 



with — 1 < A < 1, the heat transport turns out to be of ballistic type and therefore. Kg 
is infinite for infinite-length systems. Anomalous behavior of energy transport in Heisenberg 
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S" = 1/2 AFM XXZ chains has been predicted a long time ago |3i4| . It was recently associated 
with the integrability of the underlying model, because the energy density appears to be one of 
the numerous constants of motion of integrable systems Very recently, thermal transport 
in integrable models of several ID spin systems was analyzed employing the Kubo formalism 
and adopting the notion of a "thermal" Drude weight with a similar meaning as the Drude 
weight in the standard theory of electrical transport jHICI ■ The thermal conductivity can thus 
be written as |7| 

«;.(T) = «*^r)T, (2) 

where k*'* is the thermal Drude weight (in ref. |Q the thermal Drude weight is denoted as 
k = TTK*^) and r is the lifetime of the eigenstates of the spin system. The Drude- weight 
can be calculated exactly, e.g., by the Bethe Ansatz method. For a system described by 
the Hamiltonian in eq. Q), different methods of calculation give consistent values of k^^(T) 
piHIITir^l lTnilTTj . The relaxation time r is infinite for an ideal integrable ID system but turns 
finite if the influence of external perturbations, such as defects, phonons, interchain coupling, 
cannot be neglected. Theoretical work concerning the features of r for different types of 
perturbation is scarce. Below, we describe an experimental attempt to investigate details of 
energy transport relaxation in ID Heisenberg 5* = 1/2 XXZ model materials. 

We have measured the thermal conductivity k of the S* = 1/2 spin-chain compound 
BaCu2Si207 along different crystallographic directions and at temperatures between 0.5 and 
300 K. BaCu2Si207 has an orthorhombic crystal structure with lattice parameters a = 6.862 A, 
b = 13.178 A, and c = 6.897 A ^^l- Important structural elements are chains of corner-sharing 
Cu04 tetrahedra, running along the c-axis. Cu^+ ions with spins S* = 1/2 are connected along 
the chain direction via 124° Cu-O-Cu bonds. The magnetic properties of individual chains 
can be well described by the isotropic form of the Hamiltonian of eq. (Q), i.e., A = 1 and 
J = 24.1 meV ^^l- The interchain interactions are two orders of magnitude smaller ^3] and 
lead to a 3D AFM ordering transition at T/v = 9.2 K. Below Tjv, the spins are oriented along 
the c-axis. Applying a magnetic field along the easy axis leads to an exotic two-stage spin-flop 
transition ^HI; the origin of this anomalous behaviour is not completely understood |16lll7j . 
Recent interest in BaCu2Si207 was also stimulated by the speculation that below Tjy the 
"longitudinal mode" , a long-lived spin wave with a polarization parallel to the direction of 
the ordered moments, might be observed [TS|. 

Three bar-shaped samples with typical dimensions of about 0.4 x 0.5 x 2mm'^ were cut 
from the same single crystal grown using the floating-zone image-furnace technique. For each 
sample, the longest dimension was oriented along one of the main crystallographic axes. The 
thermal conductivity was measured employing the standard uniaxial heat flow technique. 
The heat input into the sample was provided by the Joule heat of a Ru02 resistor. The 
temperature gradient was monitored using a system of Chromel-Au-f 0.07% Fe thermocouples 
in the temperature region between 2 and 300 K, and by a couple of Ru02 thermometers in 
the temperature regime between 0.5 and 3 K. 

The k(T) data measured along the a, b and c-axes are shown in Fig.^ Between 2 and 32 K, 
they are in qualitative agreement with unpublished results of Takeya et al. |19| but our k- 
values are somewhat higher and the distinct features around = 9.2 K, clearly visible in the 
insets of fig. ^ are absent in the data of ref. ,19j . We suspect that this is due to the improved 
quality of the samples that were used in our study. Dips in k{T) upon magnetic ordering 
have often been observed in magnetic materials, and can be attributed to the scattering of 
phonons by enhanced fluctuations of the order parameter |ll2j . Below about 2 K, the thermal 
conductivity is almost isotropic and the temperature dependence is k oc T^'^. The anisotropy 
of k{T) is also weak above about 170 K where, for all directions, k oc T^*'-^. The anisotropy 
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Fig. 1 - Temperature dependences of the thermal conductivities of BaCu2Si207 along the a, h, and 
c axes. The insets emphasize the «:(T) features near the 3D ordering transition Tjv = 9.2 K. 

ratio / k"" for heat transport perpendicular to the chain direction is close to 1 and almost 
temperature-independent across the entire covered temperature range. It seems reasonable to 
assume that perpendicular to the spin chain direction, the phonon heat transport dominates 
and therefore, the weak anisotropy can easily be attributed to a weak anisotropy of the elastic 
constants and uncertainties in the evaluation of the sample geometry. 

At intermediate temperatures the thermal conductivity Kf^ along the chain direction signifi- 
cantly exceeds those along the perpendicular directions. At 20 K, nf^ / k"" = 2.3 and k^/k^ = 2.1. 
This particular and temperature-dependent anisotropy cannot be attributed to anisotropic 
phonon heat transport. Considering the layered crystal structure of BaCu2Si207 with CUO4 
layers separated by Ba and Si04 layers, which are stacked along the 6-axis, it is expected that 
the phonon heat transport is weakest along the ^-direction. Since this is not reflected in our 
observations, the anisotropy of the phonon heat transport is obviously very weak. A more 
detailed discussion of this point may be found in ref. ^| . 

The excess contribution to n along the c direction is most naturally associated with the ID 
energy transport by itinerant spin excitations (spinous). In electrically insulating materials, 
the total measured thermal conductivity along an a-direction is K"(r) — -|- k", where Uph 
and Ks are phonon and spin contributions, respectively. For directions perpendicular to the 
direction of weakly interacting spin chains, k", ~ because of the vanishingly small spinon 
velocities. To analyze the temperature dependence of (simply denoted as Kg below), we 
assumed that the phonon contribution Kph is nearly isotropic and subtracted Kph ~ {k°' + k^)/2 
from the total measured n'^{T). The resulting values of Ks{T) are shown in fig. [3 Since the 
measured values of k have maximum errors of about ±10% along all directions, the related 
absolute maximum of possible uncertainties in Ks(T), which are indicated by the shaded area 
in fig. El are very large at high and low temperatures, where the (T) curves for directions 
parallel and perpendicular to the chains differ only little. That is why we did not analyse 
the thermal conductivity in the AFM ordered state and restrict our analysis of Ks{T) to the 
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limited temperature region between T/v and 80 K (see dashed vertical lines in fig. 121 ■ 

Another reason for the observed anisotropy of k might be an anisotropy in the scattering 
strength between phonons and spinons. Although such a scenario cannot be completely ruled 
out, two facts suggest that it is very unlikely Firstly, the nearly isotropic heat transport at the 
highest temperatures reached in this study is not compatible with a strong anisotropy of the 
phonon thermal conductivity and thus, the phonon-spinon interaction. Secondly, the observed 
reductions of the thermal conductivities at T/v, directly related to the spin-lattice interaction 
PP, are approximately the same for all crystallographic directions. Consequently, the features 
of k(T) at T/v, clearly seen in the insets of fig. ^ are removed by subtracting + k'')/2 
from nf^ . All this suggests that, although the scattering of phonons by spin excitations may 
be relatively effective, Kph remains, nonetheless, almost isotropic and hence the anisotropy of 
the measured k(T) is almost certainly due to heat transport by spin excitations along the c 
direction. 

From the values of Ks(t) a-nd using eq. we calculated the mean free path of the spin 
excitations ts = VsT, where Vs{T) denotes the temperature-dependent average velocity of 
propagating spin excitations. Assuming that the excitations of the Hamiltonian in eq. 
are spinons obeying Fermi-Dirac statistics and adopting an energy dispersion |21| e{k) = 
( J7r/2)sin /cc, where c = 3.43 x 10^^" m is the distance between neighboring spins in a chain, 
Vs{T) can be calculated as 



de 
dk 



f{e,T)dk 



I 



fie,T)dk 



(3) 



where f{e,T) = {exp{e/kBT) + 1)"-^. The temperature dependence of the thermal Drude 
weight D*^{T), taken from ref. [HI, served to calculate r(T). The resulting values of is{T) are 
shown in Fig.|21 The shaded area again denotes the uncertainty. We note a rapidly increasing 
(■s{T) with decreasing temperature and a trend to saturation at low temperatures. 

It is rather instructive to compare this ^s(T) with available data for other similar insulating 
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Fig. 3 - Mean free path of spin excitations in BaCu2Si207 in comparison with corresponding data 
for the spin-chain compounds Sr2Cu03 and SrCu02 |20| . and the two-leg spin-ladder material 
CagLa5Cu2404i 23- The intrachain spin-spin correlation length 5, shown as the broken line, is 
from ref. |2b| 



compounds. Apart from BaCu2Si207, mean free paths of spin excitations have been evalu- 
ated for the S = 1/2 spin-chain compounds SrzCuOs and SrCuOj EHI, CuGeOa |121|2S1E1I, 
and also for the undopcd spin ladder compound Ca9La5Cu2404i. The data for Sr2Cu03 and 
SrCu02 (data from .20,), and Ca9La5Cu2404i |^ are shown in Fig.|21 For CuGeOs, the tem- 
perature dependence of £s was not calculated, but it was noted that £s does not exceed the 
value of 1300 A |23I24| . It is remarkable that aU is{T) curves in fig.Oreveal a similar trend, in 
particular the low-temperature saturation value is always of the order of 10"^ A. Also the high 
temperature features are similar and reveal a power law type T~" behavior with n > 1. It is 
important that, in spite of very different excitation spectra and therefore different tempera- 
ture variations of D*^{T) for spin chains and spin ladders (HIEI, the corresponding extrinsic 
parameters ^s(T) exhibit similar features. This suggests that the relaxation mechanisms for 
spin excitations are the same in all these compounds. 

In view of the present lack of theoretical descriptions of relaxation processes in spin chain 
compounds, only qualitative suggestions can be made. The first straightforward result is that 
the mean free path of spinons is not limited by the spin-spin correlation length. The intrachain 
correlation length in spin chains varies as ^ oc with small logarithmic corrections |2f)j . 
The result of calculations for BaCu2Si207 is shown as the dashed line in Fig.|31 The absolute 
values of £s at low temperatures are considerably larger than ^, and we also note a much faster 
temperature-induced reduction of isiT) than of ^(T) at elevated temperatures. 

In ref. ^7\, a phenomenological approximation for the spinon mean free path was postu- 
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Table I - Parameters of the fitting of Is (T) data to eg. {TP 



Parameter 
Lsd (nm) 
Asp (10^ m-i) 
T* (K) 



SraCuOs 
81 ±2 
7.1 ± 2 
177 ±5 



SrCu02 
242 ± 2 
6.7 ± 1 
204 ±2 



Ca9La5Cu2404i BaCu2Si207 



325 ±1 216 ±2 

3.3 ±0.2 41.4 ±2 

532 ± 10 80 ± 2 



lated as 



T*/T) ± L 



-1 

' sd ' 



(4) 



The first term on the right-hand side reflects the scattering of spinons by phonons, with 
the parameter Asp characterizing the spin-lattice interaction strength and T* representing a 
threshold energy typical for Umklapp-processes. The second and constant term is associated 
with the influence of magnetic defects due to non S — 1/2 ions on the Cu sites, with Lgd as 
the mean distance between them. Eq. Q is equally well applicable to all data sets shown in 
Fig. El The solid lines in Fig. O are fits to eq. with the parameter values given in Tabled 
If Lsd indeed reflects the limitation of the mean free path by defects, the fit values of this 
parameter given in Table are the consequence of similar densities of these defects in these 
materials. The exception is Sr2Cu03 where the smaller value of Lg^ might be due to the 
additional influence of bond-defects 57 ■ 

The relatively large value of the parameter Asp for BaCu2Si207 indicates enhanced spinon- 
phonon scattering in this compound. This can be understood by taking into account that for 
this material, the intrachain exchange constant J/fcs = 279 K is close to the value of the 
Debye temperature 9d which, from low-temperature specific heat data presented in ref. 28"', 
is approximately 290 K. The near identity of the energy scales for spin and lattice excitations 
leads to stronger scattering processes between the two types of quasiparticles. For the other 
compounds, the values of J/fcs are about an order of magnitude large than 9d. 

If the strongly increasing £s (T) with decreasing T is indeed due to the reduction of spinon- 
phonon scattering, it may happen that at low enough temperatures the heat input into a 
sample will produce a temperature gradient in the crystal lattice only. Correspondingly, the 
spin system adopts a constant average temperature and does not participate in the energy 
transport processes '2^. In ref. OUI we estimated that for Sr2Cu03 and SrCu02 this happens 
quite abruptly in the temperature region between about 10 and 20 K. However, the wavevector 
dependence of the spinon-phonon interaction was not taken into account. Such a dependence 
would lead to a broadening of the transition from the high-temperature region, where Hs can 
be measured by our methods, to the low-tcmpcrature regime, where Hs in not observable in 
this way. In principle, the low-temperature saturation of is{T) may reflect such a transition, 
instead of the influence of defects. 

Although cq. Q captures the main features of the temperature variation £s{T), it is not 
properly justified theoretically. A rigorous quantum mechanical analysis of the infiuence of 
the spin-lattice interaction on the energy relaxation in quantum spin chains would definitely 
be of great help. 

In conclusion, the analysis of the experimental results of anisotropic heat transport in 
BaCu2Si207 allowed for the evaluation of the mean free path ts{T) of itinerant spin excita- 
tions. A comparison with several other compounds containing isotropic Heisenberg 5 = 1/2 
chain-type structural elements demonstrates that ts (T) of all these materials exhibit similar 
features, and indicates that a universal type of energy relaxation for spinons is dominant in 
these materials. 
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